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ABSTRACT 



Context. The high-mass microquasar Cyg X-l, the best-established candidate for a stellar- mass black hole in the Galaxy, 

has been detected in a flaring state at very high energies (VHE), E > 200 GeV, by the Atmospheric Cherenkov Telescope 

MAGIC. The flare occurred at orbital phase (j> — 0.91, where <f> — 1 is the configuration with the black hole behind the 

companion high-mass star, when the absorption of gamma-ray photons by photon-photon annihilation with the stellar 

field is expected to be highest. 

Aims. We aim to set up a model for the high-energy emission and absorption in Cyg X-l that can explain the nature 

of the observed gamma-ray flare. 

Methods. We study the gamma-ray opacity due to pair creation along the whole orbit, and for different locations of the 

emitter. Then we consider a possible mechanism for the production of the VHE emission. 

Results. We present detailed calculations of the gamma-ray opacity and infer from these calculations the distance from 

the black hole where the emitting region was located. We suggest that the flare was the result of a jet-clump interaction 

where the decay products of inelastic p — p collisions dominate the VHE outcome. 

Conclusions. We are able to reproduce the spectrum of Cyg X-l during the observed flare under reasonable assumptions. 

The flare may be the first event of jet-cloud interaction ever detected at such high energies. 

Key words. X-rays: binaries - gamma-rays: theory - radiation mechanisms: non-thermal - stars: winds, outflows 
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1. Introduction 

Five X-ray binaries have been detected in the very high- 
energy region of the electromagnetic spectrum, E ~ TeV. 
Three of them, PSR B1259-63, LS I +61 303 and LS 5039, 
have been detected at different orbital phases and show 
variable emission. Four gamma-ray flares were detected by 
the AGILE satellite from the exceptional X-ray binary Cyg 
X-3 (Tavani et al. 2009). The Fermi Large Area Telescope 
(LAT) has also detected a variable high-energy source coin- 
ciding with the position of Cyg X-3 (Abdo et al. 2009). The 
fifth source, Cyg X-l, has been detected only once during a 
flare episode. This latter detection constitutes the first evi- 
dence of very high-energy gamma-ray emission produced in 
the surroundings of a stellar-mass black hole (BH) in our 
galaxy (for further discussion see Paredes 2008). 

Recently, Albert et al. (2007) reported the results from 
observations of Cyg X-l at very high energies, E > 200 
GeV, performed with the Major Atmospheric Gamma 
Imaging Cherenkov (MAGIC) telescope. No persistent 
emission was detected, but a fast transient episode was. 
The satellites INTEGRAL and Swift/BAT detected with 
some delay a related flare at hard X-rays, while only a 
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statistically poor detection was found in the RXTE/ASM 
data at soft X-rays. This wavelength-dependent behavior 
may suggest that different emitting regions were involved. 
The gamma-ray excess occurred at orbital phase 4> = 0.91. 
This can help to set constraints on the location of the emis- 
sion region. More recently, the flaring nature of Cyg X-l in 
gamma rays has been confirmed with the AGILE satellite 
(Sabatini et al. 2010). This work is devoted to a study of 
the absorption of high-energy photons in Cyg X-l and the 
implications of the resulting constraints. 

The paper is organized as follows: in the next section we 
describe the main characteristics of the source under study. 
Section 3 deals with the gamma-ray opacity by pair creation 
in the stellar radiation field. The production mechanism of 
the flare emission is then examined in the context of exist- 
ing models (e.g. Bosch-Ramon, Romero, & Paredes 2006; 
Romero et al. 2003). In particular, we explore the physical 
conditions required by the energy budget and spectrum of 
the flare event. In Sect. 4 we present a simple modelization 
for the non-thermal emission and compare our calculations 
with the observational results. Finally, in Sect. 5, we present 
a brief discussion and the conclusions. 



2. Cygnus X-l 

The binary system Cyg X-l is composed by a massive star 
and a compact object. The X-ray and radio monitoring of 
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the source over the last decades have shown that Cyg X- 
1 is most of the time in a hard X-ray state and powers 
collimated jets (e.g. Stirling et al. 2001), which makes it 
a confirmed high-mass microquasar (HMMQ, Mirabel & 
Rodriguez 1999). It is located at a distance of 2.2±0.2 kpc 
(Ziolkowski 2005). The massive star is an 09.7 lab of 40±10 
Mq and the compact object is the best-established candi- 
date for a stellar-mass BH in the Galaxy, with 21±8 M Q 
(Ziolkowski 2005). The orbit of the system is circular, with 
a period of 5.6 days and an inclination between 25° and 
65° (Gies & Bolton 1986). At radio wavelengths, a semi- 
ring surrounds Cyg X-l. This feature is thought to be the 
result of a strong shock at the location where the jet im- 
pacts onto the ambient interstellar medium (Gallo et al. 
2005). 

Regarding the flare event at VHEs, the observed energy 
spectrum is well fitted by a relatively soft power law (Albert 
et al. 2007) 



dA 



AA&t&E 



(2.3 ±0.6) x 10 _12 ( 



cm-V^eV -1 . 



E 



-3.2±0.6 



ITcV' 



(1) 



The star provides an intense radiation field that can 
absorb gamma-rays by pair creation within the binary sys- 
tem. The detection by MAGIC occurred near the superior 
conjunction, when this opacity to gamma-ray propagation 
from a region close to the compact object is expected to be 
maximum. 

The massive star has a strong wind. Considerable obser- 
vational evidence supports the idea that winds of high-mass 
stars are clumpy (e.g. Owocki & Cohen 2006, Moffat 2008). 
In a HMMQ, some clumps could eventually penetrate into 
the jet of the system enhancing the non-thermal emission, 
as proposed by Owocki et al. (2009). 



3. Gamma-ray opacity due to e H 
the stellar radiation field 

3.1. Calculations 



pair creation in 



In a HMMQ the radiation field of the massive star pro- 
vides soft photons that can annihilate gamma-rays by pair 
creation: 7 + 7 — > e + + e~ . We consider the opacity treat- 
ment for gamma-ray absorption in a massive X-ray binary 
system as in Dubus (2006) and Romero et al. (2007). The 
differential opacity for a gamma-ray at P traveling in the 
direction given by e 7 due to photons of an energy e emitted 
at S in the direction e* is (Fig. Q]) 



dr 77 = (1 — e 7 e*)n e (7 77 dedS!dZ, 



(2) 



where dfi is the solid angle of the surface that emits the 
photons and n e is the specific radiation density. 

The cross-section for photon annihilation is (Gould & 
Schreder 1967) 



er 77 (/3) 



-(I-/? 2 ) 
2(3((3 2 



>) + (3-/? 4 )ln(i±|) 



(3) 
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where /3 = (1 — 1/s) 1 / 2 , and s = E 1 e(l — e 7 e^)/(m e c 2 ) 2 
Here, E 1 and e are the energies of the gamma-ray and the 



stellar photon, respectively. This reaction occurs above a 
minimun energy given by 



2(m e c 2 ) 2 
(1 - e 7 e*) 



(4) 



where e 7 is a unitary vector in the direction of the gamma- 
ray propagation and e* is also a unitary vector in the direc- 
tion of the stellar photon propagation. The optical depth 
is a trajectory integral for which the angular dependence 
has a very significant effect. The absorption is then highly 
modulated by the orbital motion. It depends also on the 
target photon field, which is strongly anisotropic along the 
gamma-ray path . 




Fig. 1. Gamma-ray photon at P travels in the direction 
given by e 7 . This photon can be absorbed by photons of 
an energy e emitted at S in the direction e*. Adapted from 
Dubus (2006). 



Because the massive star completely dominates the 
spectral distribution of the radiative field at low energies, 
any other source of radiation for the production of pairs 
with gamma rays is neglected here. The star has a radius 
i?*, and for simplicity we asume a blackbody density radi- 
ation of a temperature T+: 



2e 2 



1 



h 3 c 3 (exp(e/fcT*) - 1) 



phem erg sr 



(5) 



The geometry considered for the gamma-ray absorption 
is shown in Fig. [3] If emission occurs at a height h above the 
compact object and perpendicular to the orbital plane, the 
distance d from the star becomes d = ydg + h? and the 
initial angle changes from ipQ = sin(0) sin(z) to ipi. Note 
that according to Fig. [2] 



cos?/>i = e r e obs 



1 



^(dl + h?) 



(do cos 2ir<fi sin i—h cos i) . (6) 
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Fig. 2. Sketch of the geometry considered for the gamma- 
ray absorption of a photon that is produced above the com- 
pact object. 

Table 1. Model parameters. Those related to the absorp- 
tion are listed first. 



Parameter [units] 


values 


R* 


Stellar radius [cm] 


1.5xl0 12 


T, 


Stellar temperature [K] 


3 x 10 4 


'"orb 


Orbital radius [cm] 


3.4 x 10 12 


e 


Viewing angle 


tt/6 


AfBH 


BH mass [M Q ] 


20 


ho 


Jet initial point [Rsch] 


50 


hint 


Height above compact object [cm] 


10 13 


Q 


Equipartition parameter 


0.1 


Fjct 


Jet bulk Lorentz factor 


1.4 


B 


Magnetic field [G] 


0.9 


V 


Acceleration efficiency 


0.1 


Ljet 


Jet kinetic power [erg s _1 ] 


10 37 


a 


Hadron-to-lepton energy ratio 


100 


qrel 


Jet content of relativistic particles 


5% 


Rjct 


Jet radius [hint] 


0.1 


e 


Thickness of the "one zone" [hint] 


0.05 


c 


Particle injection index 


2.8 


M* 


Stellar mass loss rate [M yr -:L ] 


3xl0" 6 


Vex, 


Terminal wind velocity [cm s _1 ] 


2 x 10 8 



The parameters adopted for the calculations are shown 
in Tabled! 

Under adequate conditions, the absorption, resulting in 
the creation of energetic pairs, and the Inverse Compton 
(IC) emission from them, can operate in an effective way 
to develop electromagnetic cascades which can consider- 
ably modify the original gamma-ray spectrum (see e.g. 
Bednarek 1997 and Orellana et al. 2007 for detail treat- 
ments). Electrons with TeV energies in the stellar radiation 
field may also lead to this situation. At TeV energies the 
rate of electron energy losses in the Klein-Nishina regime 
is reduced by the diminution of the IC cross-section. The 
ambient magnetic field must be smaller than a critical value 
B c for the synchrotron losses not to overcome the IC ones. 
In order to determine if effective electromagnetic cascading 
can occur within the system it is then necessary to know 
the magnetic field strength in the gamma-ray propagation 
region. Such a field is dominated by the stellar magnetic 
field. Magnetic fields measured in massive stars can reach 
~ 10 3 G, which is much greater than the critical value B c . 
For close binaries like Cyg X-l we can expect that B > B c 
(Bosch-Ramon, Khangulyan, & Aharonian 2008) over the 
whole region of gamma-ray production. We here assume 



that B > B c , and neglect the effects of electromagnetic cas- 
cades, as well as the reprocessing of the absorbed energy by 
synchrotron radiation. The latter situation was considered 
by Bosch-Ramon et al. (2008), who deal with the diffusion 
of secondary pairs into the system. Zdziarski et al. (2008), 
on the other hand, do consider that the HE photons iniciate 
a spatiallty extended pair cascade, but we will comment on 
this below (Sect. 5). 

3.2. Results 

In Fig. [3] we show a 2D-map of the attenuation coefficient 
e~ T as a function of the energy E and the height h above 
the orbital plane. This absorption map corresponds to the 
orbital phase = 0.91, when the flare occurred. As can 
be seen from the figure, the attenuation is high at energies 
between 10 GeV and 10 TeV, close to the compact object, 
which makes the absorption problem in the energy range 
where MAGIC detected the flare very relevant. 

In Fig. [4] we show a 2D-map of the attenuation coeffi- 
cient for E = 1 TeV as a function of the orbital phase </> 
and the height h. It can be seen that the absorption drops 
strongly as the height above the compact object increases 
for h > 10 11 cm. When h < 10 11 cm the absorption does 
not present major changes, due to the distances involved 
that make the photon density remain rather constant (i.e. 
R* ~ 10 12 cm and r orb = 3.4 x 10 12 cm; see Fig. |2J). Bosch- 
Ramon et al. (2008) find out from opacity calculations near 
the superior conjunction that the TeV emitter in Cyg X-l 
should be located at a distance greater than 10 12 cm above 
the compact object. Our absorption calculations agree with 
this result. Notice that our results cover a much larger pa- 
rameter space. 




Log E [eV] 

Fig. 3. Absorption map as a function of the height h above 
the compact object and the energy E for orbital phase <\> = 
0.91. 



From Albert et al. (2007) the observed flux is a power 
law (Eq. [TJ, in the energy range between 150 GeV and 3 
TeV. Considering that the intrinsic flux from the flare is also 
a power law Fi nt — nE~ a , we can relate both expressions 
through 



F, 



obs 



-Flute 



-t{E) 



(7) 
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Fig. 4. Absorption map as a function of the orbital phase 
4> and height h above the compact object for energy E = 1 
TeV. 

From the computed numerical values of t(E), using the de- 
pendence of r on the height h, we obtain the intrinsic spec- 
tral index a as a function of the latter parameter. Figure [5] 
shows the result. Note that for an altitude of h ~ 10 13 cm 



6 




Fig. 5. Range within the error bars of the intrinsic flux 
index as a function of the height h. 



the de-absorbed and production spectra are essentially the 
same. 

4. Flare production mechanism 

A hadronic MQ model for Cyg X-l has been already con- 
sidered in Orellana et al. (2007) based on ideas advanced 
by Romero et al. (2003). We here revisit that scenario with 
the addition of the interaction between the steady jet and 
a more dense target: a clump from the stellar wind that 
allows through locally generated shocks the reacceleration 
of the particles that produce VHE emission far from the 
BH, as in Araudo, Bosch- Ramon, & Romero (2009). The 
jet+clump system is assumed to be momentarily in steady 
state. As observed in the stable configuration of a micro- 
quasar in a low-hard X-ray state (e.g. Fender, Belloni, & 



Gallo 2004), we assume a continuous jet. The calculations 
of the emission are based on the works by Bosch- Ramon et 
al. (2006) and Romero & Vila (2008). 

The jet is considered perpendicular to the orbital plane, 
and launched at a distance ho above the compact object. 

We consider that farther down the jet the magnetic field 
reaches values well below equipartition. Following Bosch- 
Ramon et al. (2006) the magnetic field in the jet reference 
frame can be calculated as 

B(h) = v/^8^;. (8) 

In Eq. (JSJ) g> is the equipartition parameter and e p is the 



matter energy density. Then, 



'jet 



irRt 



jet' l; Jet TO P 



-(-Epl 



'jot 



2nh 2 



Vjet) 



(9) 



where Wj et is the bulk velocity of the outflow, we set Uj ct ~ 
0.7c (Heinz 2006). The jet radius is i?j e t = xh, and (jE7 p k) 
is the mean kinetic energy of the cold proton, taken to 
be the classical kinetic energy with a velocity equal to the 
expansion velocity of the jet (v CKp = x w jot)- 

A small fraction of the jet power is transformed into 
rclativistic particles in a "one-zone" acceleration region lo- 
cated above the compact object, at the height of the impact 
with the clump. Here we assumed h- lnt = 10 13 cm, based on 
our opacity constraints. 

The kinetic power in the form of relativistic particles 
is assumed to be proportional to the jet's power, L xc \ = 
<7 re iLj e t, with q IC \ = 0.05 and Lj et = 10 37 erg s _1 (Gallo et 
al. 2005). We considered both hadronic and leptonic con- 
tent, L rc \ = L p + L c . The ratio of relativistic protons to 
electrons luminosity a in the jet is unknown. We adopted 
a = 100, a similar value to what is observed in the galactic 
cosmic ray spectrum (e.g. Berezinskii et al. 1990). 

The minimum kinetic energy is taken to be on the order 
of the rest mass energy of the corresponding particles. The 
maximum energy for the electrons is obtained equating the 
cooling rates with the acceleration rate. The acceleration 
rate by Fermi mechanism, £~ c * = E~ 1 dE/dt, of a particle 
with energy E in a magnetic field B, is given by 



t" 1 = 

ace 



r/ecB 



(10) 



with rj the acceleration efficiency, which is assumed here 
to be high, ~ 10 %. The maximum energy for protons is 
restricted by the size of the acceleration region because 
the particle giroradius r g = E/eB should not exceed i?jet- 
The energy losses considered for electrons are adiabatic, 
IC, synchrotron and relativistic Bremsstrahlung, and are 
calculated in the jet reference frame (RF). 



For adiabatic losses, the cooling rate is 



'ad 



__ 2 fjet 

3 ftdnt 

The synchrotron losses rate is 



-l 



4 (jtcUb fin 



synchr 



E 



3 rripC 2 \ m J mc 2 ' 



(11) 



(12) 



where or is the Thomson cross-section and Ub is the mag- 
netic energy density. 
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The IC loss rate can be calculated from (Blumenthal 
Gould 1970) 



r _1 — — 



' IC 



E„ 



hE v 



( A dN A 

( ei ~ e >17T~ dei ' 



(13) 



where e and ei are the incident and scattered photon ener- 
gies, respectively, and 



dN _ 1 27rr§mc 3 n ph (e)de 



didei E c 



with 



/(<?), 



f(q) = 2qlnq+(l + 2q)(l-q) + 



1 (bq) 2 



2a + bq 



(i-«). 



(14) 



(15) 



Here b = Ae'y/mc 2 and q = ei/[b(E e — ei)]. The seed photon 
field is provided by the companion star, considered as a 
blackbody at T* and is transformed to the jet reference 
frame (Dermer & Schlickeiser 2002) 



n e o 



I? ct (l + /3M) 2 



(16) 



where Fj et is the jet Lorentz factor, Q represents the photon 
direction, \x — cos 0, and O is the angle between the photon 
direction and the jet axis (the quantities with primes are in 
the jet RF). We considered the "head on" approximation, 
in which fj, = — 1. 

The relativistic Bremsstrahlung losses for a complete 
ionized plasma were computed according to (Berezinskii et 
al. 1990) 



Br 



4nZ r„ac 



In- 



2E f . 



1 



m G cr 



(17) 



where n = n w ind/r is the density of target ions expressed 
in the jet RF. For external Bremsstrahlung the target-ion 
density is that of the stellar wind ions (the wind is consid- 
ered as a completely ionized plasma). At a height h 



n w [ n d(h) 



M* 



4:TTv oc ,m p (h 2 + a 2 ) 



R* 



VTP 



r rb z 



, (18) 



where Voo is the terminal velocity of the wind (e.g. Romero 
et al. 2003 and references therein). In order to take into ac- 
count the mixing between the jet and the background wind 
material in a phenomenological way, we introduced a pene- 
tration factor / p = 0.3 (Romero, Christiansen, & Orellana 
2005). The clump is considered as a condensation of the 
wind with a density of ~ 10 14 cm" 3 . Notice that particle 
rejection at the jet-wind boundary is considered only for 
the background wind. The clump is assumed to fully pene- 
trate the jet (see Araudo et al. 2009 for details). For internal 
Bremsstrahlung the target-ion density is the proton density 
of the jet, directly derived from rhj et = Lj et /Tc 2 . 

Relativistic protons lose energy through adiabatic ex- 
pansion, synchrotron radiation, and by losses produced by 
hadronic interactions. The energy loss rate produced by 
proton-proton interactions is 



pp ^p ^^pp pp ' 



(19) 



where n p is the density of target protons and K pp the in- 
elasticity (~ 0.5). The cross-section can be approximated 
(Kelner, Aharonian & Bugayov 2006) by 



cr pp = (34.3+1.88X + 0.25L 2 ) 



i.h 



E,„ 



vab. 



(20) 



where L = ln(£ , p /lTeV). Photomeson production is not 
considered because the stellar photons do not have enough 
energy to reach the threshold energy of this process. 
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Fig. 6. Acceleration and cooling rates at /i; n t = 10 13 cm in 
the jet for primary electrons and protons. 

In Fig. [5] we show the rates of cooling and energy gain 
for electrons and protons in the acceleration region, which 
is considered to be the bow-shock between the jet and the 
clump. The electrons reach TeV energies while the protons 
can attain energies ~ 10 2 TeV. 

In the one zone approximation the steady state parti- 
cle distributions N(E) result from the transport equation 
(Ginzburg & Syrovatskii 1964) 



d_ 

~dE 



&E 



At 

where t, 



N(E) 



loss 



N(E) 



to 



Q(E), 



(21) 



hint /V jet- 

The exact analytical solution of the equation is 

l 



N(E) 






r-E 



/* AE' Q(E') 



loss 
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xcxp 



t(E,E>) 



with 
t(E,E>) 



E' 



dE' 1 



dE" 

~dT 



(22) 



(23) 



less 



The particle injection function, Q(E), is assumed to be a 
power-law in the energy of the particles, 



with 

F{q) =2q\n q +{l + 2q){\ - q) + \{l - l)j^^-y (30) 

and n = (4eE c )/{(m c c 2 ) 2 ) , q = {E 7 )/[(QE e (l - E n /E e ))]. 
The relativistic Bremsstrahlung contribution is given by 



Q(E) = Q E 



-C 



(24) 



This distribution is expected to be the result of diffusive 
particle acceleration by the reverse shock. The index ( is 



where 



E n V 



na B (E e ,E^—N e (E e )dE e , 



Air 



Aar 2 . 



taken as 2.8 for both types of particles, in accordance with (T b( 1 -c<c, -by) — —j^ <P\^c, -by), 
the obtained results shown in Fig. [5] The normalization 
constant Qq for each type of particle is derived from L CtP 
as 



L e , p = V 



^-^G^p^Cjp^CGjpx^Gjp)-! 



(25) 



E~ 

and 

cf>(E e ,E^= [1 + (1-E y /E c ) 2 - 
2E C (E C — Ej 
m c c 2 E^ 



(31) 



(32) 



x On 



2/3(1 - E 7 /E e )] 
~ 2 



(33) 



All luminosities were calculated in the jet co-moving 
RF. Photon energies in both frames are related by the 
Doppler factor D as 



(34) 



where V is the co-moving one-zone volume. 

4.1. Radiative processes 

We consider synchrotron emission from both electrons E~ ( = DE^ 
and protons, inverse Compton emission from electron 
interactions with the stellar photon field, internal and where 
external relativistic Bremsstrahlung, and inelastic colli- 
sions between relativistic protons in the jet and the cold D 
material that forms the jet, plus with the matter of the 
clump and the background wind. We checked that the The i uminosity in the observer frame is given by (e.g. Lind 
emission produced by secondary particles is negligible, as ^ Blandford 1985) 
well as the synchrotron self-Compton (SSC). 



r jet (l-/3j et cos6'p)' 



(35) 



The synchrotron emission was computed with the ap- 
proximation 

^(Ey)= E,V^i- /*"" dEN(E)^lM 



L 7 (£ 7 ) = D 2 L'JE') 



(36) 
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In order to compute the gamma-ray emission produced 
by neutral pion decay we note that the p — p cross-section 
parametrization ( Eq. (|2"0|) ) is given in the laboratory RF. 
Then, we convert the flux of relativistic protons to the lab- 
oratory frame: 



where 



E, 



3 ehB ( E 



(26) 



(27) 



AT- 



J(E p ,9p) 



E p - f3 cos Op JE 2 - m 2 c A 



-in 



An mc \mc 2 
and the usual meaning for the constants c, h, e. 

The IC emission by the electron population was calcu- 
lated as 



sin(92-|-r 2 ( cos6> p 



0E P 



y/n 



1/2 



,(37) 



where A is a normalization constant. The flux of protons, 
which is isotropic in the jet RF, is beamed in the lab RF, 
as indicated by the dependence on the viewing angle 9 p . 

The gamma-ray luminosity, for E p < 0.1 TeV, can be 
obtained straightforwardly as 



Lic(E~f 



E 2 V 



dE e N e (E e ) 



dePic(£ c ,£ 7 ,e), 



JB mln y/E£ - mic 4 



(38) 



where the spectrum of scattered photons is 



v °> with E m {n — Ej + m 7r c 4 /4_E 7r . In the formalism of the 5- 
functional (Aharonian & Atoyan 2000) the 7r°-emissivity is 
given by 



P lc (E e ,E 7 ,e) 



3g T c(™ c c ) n ph ( 
4E 2 e 



-F(q), 



!1± 



(29) qir(E w ) = -^-Upp ( m p c 2 + E^/n-n ) J p I m p c 2 + E^/k^ ) (39) 
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Fig. 7. Computed SED and the MAGIC observational data from Cyg X-l (Albert et al. 2007). A two-temperatures 
corona with a non-thermal component is presented as well. The data are from McConnell et al. (2000). The similar data 
from Malzac et al. (2008) can be easily fitted (see Romero, Vieyro, & Vila 2010). 



with k* ~ 0.17 (Gaisser 1990). For E p < 0.1 TeV down 
to the threshold, a slightly modified version of ^-functional 
approximation is needed, using the replacement 



S{E K — rtTr-Ekin) — > hS(E n — K-xEkin). 



(40) 



Here h is the total number of 7r° created per p — p collision. 
The gamma-ray luminosity in the range 0.1 TeV< E p < 
10 5 TeV can be obtained from (Kelner et al. 2006) 



L 7 (£ 7 ) = n p E*V 



&inel(Ep)Jp(E p ) 



XR 



£L 



^ \ dE p 

E p ' J Ep 



.E 



with Fjl -^-^E p I a function of E 1 and E p . For further 

details on radiative processes see Vila & Aharonian (2009) 
and references therein. 

In order to reproduce the observed spectral energy dis- 
tribution (SED), the density ratio between the clump and 
the wind at the base is ~ 4.6 xlO 4 , i.e. n c ~ 3.3 x 10 14 
cm~ 3 . 

Figure [7] shows the computed SED. We have included 
the thermal emission by the star, which largely dominates 



at optical energies. At X-rays, the components of the emis- 
sion by the accretion disk and a corona should be added to 
our results. These components in the low-hard state have 
luminosities ~ 10 37 erg s _1 and extend up to ~ 150 keV (see 
Romero et al. 2002), in a way that they completely dom- 
inate over the non-thermal radiation. The emission from 
the corona and a non-thermal tail (McConnell et al. 2000, 
Malzac et al. 2008) are also shown. The model for this emis- 
sion is from Romero, Vieyro, & Vila (2010) and is presented 
in detail elsewhere. Here we show only the results relevant 
to Cyg X-l. 



(41) 4.2. Internal absorption 



Internal photon-photon annihilation within the region of 
gamma-ray production can result in strong attenuation of 
the radiation (Aharonian et al. 2008, Romero & Vila 2008). 
The opacity is again an integral of Eq. ((2|), but now con- 
sidering the locally produced photons with density n p h(e). 
We can use the symmetry in one of the angles to write 



r(E 7 ) = - 



l -'eth 



(1 - u) cr 77 (/3)n ph (e)duded/. (42) 



Here, u = cos ■&, ■& is the angle between the momenta of the 
colliding photons, I is the photon path, and the cross-section 
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<7 77 (/3) is given by Eq. Q. The absorbing photon fields are 
those generated within the jets (i.e. those calculated in the 
former section). At energies E 1 <; 10 15 eV, the dominant 
absorbing field is the synchrotron radiation from electrons. 
In the local approximation of Ghisellini et al. (1985), 



^-synchr 



^synchr 



(43) 



where e syn chr is the synchrotron power per unit volume per 
unit energy: e sync hr = L 1 /{e 2 V), with L 7 from Eq.fl^G]). 

The geometry considered requires r = i?j e t and < I < 
i?j Ct . We find that t{E 1 ) is completely negligible (at the 
level of r ~ 10~ 6 ), implying that the attenuation coefficient 
is- 1. 



5. Discussion 

The VHE transient emission of Cyg X-l occurred when 
the BH was behind the star with respect to the observer. 
Because of the high absorption in the flare detection en- 
ergy range, the emission close to the BH is not enough to 
explain the observations, unless the photons travel far away 
from the star, initiating a spatially extended pair cascade 
as considered by Zdziarski et al. (2008). This requires a 
fine tunned magnetic field, which allows the instantaneous 
isotropization of the electrons, but does not overcome their 
IC radiative losses. A more realistic/accurate calculation 
of the electromagnetic cascade propagation is then desir- 
able. Such simulations (following the electron trajectories) 
will be available in a future work as an application of the 
code developed by Pellizza et al. (2009). Previous ID cas- 
cade simulations (Orellana et al. 2007) are consistent with a 
strong absorption and steep spectrum at TeV energies. The 
results by Bosch- Ramon et al (2008) have shown that if the 
cascades are suppressed by effects of the magnetic field, the 
synchrotron emission of the secondary pairs peaks at lower 
energies (~ GeV). 

Romero, Kaufman-Bernado, & Mirabel (2002) have sug- 
gested that Cyg X-l could go through occasional microb- 
lazar phases and have estimated that the luminosity in the 
observer RF can be up to one order of magnitude higher 
than the luminosity in the jet RF. Even taking this into 
account, a flare triggered at the base of the jet is unde- 
tectable due to absorption at phase 0.91. A remaining op- 
tion could be a very short episode with a highly increased 
acretion/ejection rate, but this is speculative given the lack 
of evidence at lower energies supporting the hypothesis. 

Under the geometry considered here (a jet perpendicu- 
lar to the orbital plane, which has an inclination of — 30 
deg), the high-energy emission should have occurred at a 
large distance above the compact object where the absorb- 
ing photon field is diluted. In order to quantify the radiative 
outcome in this scenario we have considered the interaction 
of relativistic particles accelerated in a narrow region of the 
jet with the target particles of a dense clump of the wind. 

The flare timescale is related to the permanence of the 
clump inside the jet. For a spherical clump with a radius R c 
smaller than the jet radius i?j e t ~ 10 12 cm we can make a 
zerolth order estimation of the time that it takes the clump 
to cross the jet: t c . The clump velocity is the wind velocity, 
which at this height is simply u^: 



tr 



2R jct / Voo ~10 4 s. 



(44) 



The flaring episode had a timescale shorter than one day 
and a rising time of about one hour, which is on the same 
order as the t c estimated. 

The simple model presented here for the broadband 
spectrum of Cygnus X-l reproduces fairly well the observed 
SED by MAGIC during the flare using a set of param- 
eters that agrees with reasonable values for this source. 
Interactions between the clumpy winds of massive stars 
with the relativistic jets in HMMQ are expected to be pro- 
duce flaring episodes at high and very high energies, and 
may be detectable by the new high-energy detectors, like 
Fermi, MAGIC II, and VERITAS. 
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